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Abstract. The application of the Stromgren photometric luminosity calibra- 
tions to different types of CP stars is reexamined in the light of the new Hip- 
parcos data. A first attempt is made to use the LM statistical parallax method 
(Luri et al., 1996) - based on the maximum likelihood principle - to obtain a 
calibration of the absolute magnitude as a function of two Stromgren colour 
indices, thus reflecting effective temperature and evolution. Its application to 
a sample of Si stars and to a sample of normal main sequence stars in the same 
spectral range allows us to compare the calibrations obtained and to discuss 
the position of Si stars in the HR diagram. Additionally, a sample of bonafide, 
spectroscopically selected Am stars together with normal main sequence stars 
are used to derive a new absolute magnitude calibration for late A-type main 
sequence stars, taking into account the effects of evolution, metallicity and 
stellar rotation. 

Key words: Stars: chemically peculiar - Stars: fundamental parameters - 
Stars: early type - Stars: distances 

1. Introduction 

To deal with luminosity calibrations and to obtain a good exploitation of the 
high precision parallax data obtained by Hipparcos, it is required to use both, 
robust statistical methods, capable to take into account physical characteristics 
of stars (i.e. evolution, metallicity, rotational effects, etc.), and a well defined 
spectroscopic sample. Taking advantage of the fact that Stromgren photometry 
has proved to be a powerful tool for characterizing the physics of main sequence 
early type stars, some preliminary work is presented here for application to Si 
and Am stars. 

Our sample contains all CP stars included in the Hipparcos Catalogue (ESA, 
1997) having good spectroscopic information on peculiarity and complete Strom- 
gren photometry in the Hauck and Mermilliod (1996) compilation. The selection 
of the stars according to their spectral types is fully explained in Gomez et al. 
(1998, this colloquium). For comparison, a sample of normal main sequence 
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Table 1. CP stars with uvbyHp photometry. 





Selection from 
Gomez et al. (1997a) 


with uvbyHp 
Hauck & Mermilliod 
(1996) 


Range of 
SP 


Si 


440 


173 


B5-A3 


HgMn 


76 


69 


B7-A0 


Sr-Cr-Eu 


378 


172 


B8-A9 


Am 


1059 


533 


A0-A9 


Normal B0-A9 V 


3460 


1589 


B0-A9 



stars in the range B0-A9 V, selected from the Hipparcos Survey (Gomez et 
al., 1997), has been also used. Table [j] shows the stars in each CP group (the 
sample is reduced to less than 50 % when Stromgren photometry is required). 
Using the different samples, a test of the capability of the uvbyHp system to 
detect CP stars has been performed. The ([mi], [ci]) plane has been classically 
accepted as the most discriminant plane in the Stromgren system to separate CP 
from normal stars. Thus for example, Crawford (1979) established the criterion 
Smi ^ —0.020 to separate Am stars from A3-A9 normal main sequence stars. 
In contrast to Abt (1984), who concluded that about 75 % of the Am stars 
could be photometrically classified as such, we obtain that only 55 % of the 
classical Am (Sp(k)-Sp(m) ^ 5) and 33 % of the proto-Am (Sp(k) - Sp(m) < 
5) can be detected using only Stromgren colours. Philip et al (1976) criteria 
(E(b — y) ^ -0.040) to detect hot CP stars is found to be unable to acomplish 
the objective (only 2 % of Si stars and 5 % of Sr stars has been detected) . Taking 
into account that peculiar spectral features alter all Stromgren indices, Masana 
et al. (1998) established a new criterion of detection defining a Ap parameter 
- a linear combination of several uvbyHp indices obtained through Multiple 
Discriminant Analysis. Applying this criterion to the hot stars in the present 
sample (see Figure P we photometrically classify as peculiar 31 % of the Si stars 
and 56 % of hot Sr-Cr-Eu stars, whereas no HgMn star is detected. Furthermore, 
as indicated by Masana et al. (1998), the reddening decreases Ap, so reddened 
hot CP stars may be seen as non-reddened normal stars, thus decreasing its 
capability for detection. We can conclude that the Stromgren system is suitable 
but less powerful for a photometric detection of peculiarities than the Act system 
(Maitzen, 1976; specifically designed to measure characteristic features on the 
spectra of peculiar stars) or even the A(Vi — G) combination of the Geneva 
indices (Masana et al., 1998). 

A crucial point to attack the problem of the M v calibration of CP stars is 
the reddening correction. For CP2 stars, the peculiarities in the spectra result 
in bluer colour indices, the (b — y) and c\ indices being smaller, thus leading 
to underestimate their reddening when treating them as normal stars. In agree- 
ment with Adelman (1980) and Maitzen (1980), Masana et al. (1998) derived a 
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correction to the E(b — y) obtained when using standard relations valid for nor- 
mal stars as a function of the Ap parameter (AE(b —y) = —0.001 + 0.008Ap). 
Even applying this correction to the stars with Ap > 1.5, we obtain that 57 % 
of Si stars and 38 % of Sr-Cr-Eu (early region) stars have E(b — y) < , com- 
pared with only 9 % of normal main sequence stars in the same spectral range 
with negative excess (using Crawford's (1978) standard relations). The problem 
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Figure 1. Detection of hot CP stars using the Ap parameter: Ap = p — p — 1.5 
(dashed line), p = f([u — b]) being the relation for normal stars (full line) 

is maintained in the intermediate region (Grosb0l's (1978) standard relation): 
56% of Sr-Cr-Eu stars have negative excess compared to 37 % of normal main 
sequence stars. In the late region (A3-A9), where Crawford's (1979) standard 
relation is considered, the overestimation of reddening is also clear for Sr-Cr-Eu 
(18 % with E(b — y) < compared with 37% for normal stars). The same excess 
distribution is obtained for Am and normal A3-A9 stars (35 % and 37 % of 
stars with E(b — y) < respectively). Having in mind these biased photometric 
reddening determinations when peculiarities are present, we preferred to use 
the interstellar absorption model by Arenou et al. (1992) in the LM method, 
which gives reddening corrections as a function of the star position (r, /, b). Im- 
provement of the interstellar absorption models using the new Hipparcos data 
is needed. 



2. Testing the current available Mv photometric calibra- 
tions 

The photometric distances computed following the algorithms of Crawford (1978), 
Balona & Shobbrook (1984), Jakobsen (1985) and Eggen (1974) for the stars be- 
longing to the early region, are compared in Table @ with the distances derived 
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directly from Hipparcos parallaxes for Si, Hg-Mn and Sr-Cr-Eu stars as well as 
for B4-B9 V normal stars. Only stars with good parallax, IS.it/it < 30%, have 
been used. All these algorithms, based on normal main sequence stars, rely on 

Table 2. Mean differences between Hipparcos and Stromgren photometric distances 
(only stars with An/n < 30% are considered), Units: parsecs. 



Early Region 


Calibration 


Si 

(107 *) 


Hg-Mn 

(67 *) 


Sr-Cr-Eu 

hot 
(47 *) 


B4-B9 V 

normal 
(361 *) 


Crawford (1978) 

Balona & Shobbrook (1984) 

Jakobsen (1985) 

Eggen (1974) 


-16 ± 67 
-28± 84 
-15 ± 65 
-61 ± 97 


+3 ± 48 
+0 ± 50 
-4 ± 51 


+6 ± 40 
+3 ± 48 
-14 ± 55 


+ 1 ± 60 
-1 ± 60 
-5 ± 61 
-19 ± 69 


Intermediate and late Regions 


Calibration 


Sr-Cr-Eu 

intermediate 
(26 *) 


Sr-Cr-Eu 

late 
(36 *) 


A0-A2 V 

normal 
(338 *) 


A3-A9 V 

normal 
(376 *) 


Stromgren (1966) 
Guthrie (1987) 


+23 ± 43 


+62 ± 84 


+6 ± 31 


+ 10 ± 38 



the fact that the (3 index is primarily a function of luminosity and that c D and 
[u — b] are related to temperature in the early region. As the c Q index of CP2 
stars falls below those of normal stars for a given spectral type (Preston, 1974), 
systematic trends in the M v derivation can be present when using the calibra- 
tions by Crawford (1978), Balona & Shobbrook (1984) and Jakobsen (1985). 
Napiwotzki et al. (1993) and Megessier (1988) pointed out that the [u — b] index 
is slightly more sensitive to temperature than c Q and can be accepted as a good 
indicator of temperature for Si stars. Eggen (1974) used this index but his cali- 
bration shows strong systematic trends. On the other hand, the /? index is also 
affected by the peculiarities in the spectra, as pointed out by Hauck (1975). 

From Table || we see that, due to the above-mentioned effects and to the 
underestimation of reddening (when computed using Stromgren photometry), 
the photometric distances of Si stars have always been overestimated. No signifi- 
cant differences are found neither for Hg-Mn stars (as pointed out by Stromgren 
(1966) these stars behave as normal stars as far as the ci, mi and (3 indices are 
concerned) nor for B4-B9 V normal stars. The standard deviation obtained (~ 
60 pc for normal stars) agrees with a relative error in photometric distances of 
about 30%. 

Although the intermediate and late regions contain few Sr-Cr-Eu stars, Ta- 
ble shows that the photometric distances are underestimated when classical 
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calibrations for normal stars are used. Again, the overestimation of reddening 
for these stars can contribute to this effect. 

3. A new calibration for Si stars using the LM method 

The LM (Luri et al., 1996), based on the Maximum Likelihood principle, allows 
a simultaneous determination of a luminosity calibration, kinematical charac- 
teristics and spatial distribution of a given sample. In a different approach that 
the one presented by Gomez el al. (1997a) in this colloquium, we undertake a 
derivation of a calibration of absolute magnitude of Si stars and B4-B9 V normal 
main sequence stars by parametrizing M v as a linear combination of Stromgren 
colour indices (Q), the observational entries being then (a, S, m, ir, fi, V r , Q). No 
selection function for colours is imposed (S(Q) — 1), and an asymmetrical gaus- 
sian distribution of Q around a maximum (C Q ) colour index is modelled. For a 
fixed value of Q, a gaussian distribution for M v is assumed: 

1 , M — A-f q (Q) , 2 

®m(M\Q) = e~* ( «u ) (1) 

and as a first attempt we consider the dependence of M v on the effective tem- 
perature: 

M (Q) = P i+ P 2(log w ([u -b}+ 0.2)) (2) 

Pi and P2 are simultaneously determined together with the absolute magnitude 
dispersion (ctm„) , the parameters defining the velocity ellipsoid (mean (U, V, W) 
and dispersions (<Jui &v > a w)) an d the scale height (Z ). We use the [u — b] in- 
dex as explained above and a logarithmic function allowing us to solve the non- 
linearity of the main-sequence in this spectral range (Balona and Shobbrook, 
1984). The method, able to identify and characterize distinct physical groups 
in inhomogeneous samples, has been applied considering the presence of two 
groups. Results for Si stars indicate the existence of a main group (containing 
98 % of the sample) and a secondary group (2%) probably composed of misclas- 
sified objects and/or high velocity stars. Approximately the same percentage 
of classification is obtained when normal stars are considered. The relations 
obtained for the main groups are: M = —0.10 + 4.02(logi ([u — b] + 0.2)) for 
normal B4-B9 V (416 stars) and M = 0.01 + 3.58(log 10 ([u - b] + 0.2)) for Si 
stars (149 stars). 

One of the major advantages of the Stromgren system, compared with other 
photometric systems for early type stars, is its capacity to account for evolu- 
tionary effects through the (3 index. A first attempt was made to incorporate 
directly in equation (|J) a term depending on f3 but the high correlation between 
P and [u — b] did not allow to obtain stable values for the Pi coefficients. As a 
preliminary attempt we have applied the LM method in two steps: once the LM 
algorithm has been applied using equation (^|), we seek for a second solution 
considering: 

M - M (l srt step) - P [ + p' 2 {(3) (3) 
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For the main group (98 % of stars) the following relations have been obtained: 

M = -16.44 + A.02(log 10 ([u -b] + 0.2)) + 5.83/3 (4) 

for normal B4-B9 V (416 stars), with <j^i = 0.60 mag (in the l rst step : 0.70 
mag), {U .V , W ) = (-12,-15, -7) km/s, {au,a v ,a w ) = (9,11,5) km/s, Z a = 63 
pc, and 

M = -21.70 + 3.58(/o 5 i ([u - b] + 0.2)) + 7.86/3 (5) 

for Si stars (149 stars), with ou = 0.77 mag (in the l srt step : 0.83 mag), 
{U .Vo,W ) = (-15,-15, -8) km/s, (au,a v ,a w ) = (9,10,6) km/s, Z = 66 pc. 
Relations (4) and (5), plotted in figure || for different (3 values, are defined for 




[u-b] 

Figure 2. Application of the LM method in two steps in which the evolutionary effects 
are acounted for. Full lines: relations obtained for normal B4-B9 V stars plotted at (3 
= 2.6, 2.7 and 2.8 (from top to bottom of the figure). Dashed lines: the same for Si 
stars 

the whole population, so the mean relations are corrected for the Malmquist bias 
which is present in the samples used, which are limited by apparent magnitude. 
We can see a different behaviour of the (3 index as indicator of luminosity for 
normal and Si stars, the differences being more pronounced when decreasing the 
effective temperatures. Nonetheless, from the position of individual stars in the 
HR diagram computed with this solution, we can conclude that the evolutionary 
state of Si stars is nearly the same as that of normal main sequence stars. More 
work will be devoted in the near future to finding, through Principal Component 
Analysis, a suitable proper combination of colour indices that could allow a 
simultaneous determination of all the pi parameters using the LM method. 
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4. A new calibration for Am stars using the BCES method 

A different statistical treatment is used for the determination of the luminosity 
calibration in the late Stromgren's region, taking into account the dependence 
of M v on effective temperature, evolution, mctallicity and rotation. BCES fit- 
ting algorithm (Bivariate Correlated Errors and intrinsic Scatter, Akritas & 
Bcrshady, 1996) is an extension of the ordinary least squares estimator for two 
variables. It allows the consideration of measurement errors on both variables 
and the use of correlated errors between them, and also takes into account un- 
known intrinsic scatter and the possibility that the size of the measurement 
error depends on the variables. 

This algorithm has been applied to a carefully selected sample of 97 Am 
stars (North et al., 1997) and very accurate Hipparcos trigonometric parallaxes 
(Att/t: ^ 0.14). The photometric data has been taken from Hauck and Mcrmil- 
liod (1996) and the projected rotational velocity Vsini from Abt and Morrell 
(1995). Individual correction for binarity has been applied whenever possible 
(30 % of the sample) and a mean correction of 0.2 mag for SB1 binaries and 
stars with variable radial velocity has been considered, leaving 29 stars without 
correction. No correction for interstellar absorption has been applied, after re- 
alizing that it is meaningless when comparing it with the uncertainties in the 
binarity correction, and considering that the E(b — y) computed using Craw- 
ford's (1979) standard relation is distributed around a mean value of 0.0014 
mag. 

To estimate the contribution of the stellar rotation and the blanketing effects 
(which are highly correlated in this spectral range) on the photometric M v 
calibration, it is necessary to have a large range of Vsini and Smi. For this 
reason, the calibrations have been obtained by joinig to the sample of Am stars 
a sample of normal A3-A9 V stars carefully selected to be non-binaries and non- 
peculiar. The quality and source of its trigonometric parallaxes, photometry and 
Vsini is the same as for the Am stars. 

An evaluation of the existing photometric calibrations of luminosity through 
their comparison with the Hipparcos data, allows us to conclude that Crawford's 
(1979) calibration is systematically shifted by 4 mas in parallax {^Hipp < ^Craw) 
for Am stars, but shows no systematic trend for normal A3-A9 stars. The inclu- 
sion of 8m\ and Vsini terms proposed by Guthrie (1987) reduces Crawford's 
(1979) shift for Am stars but introduces a systematic deviation for A3-A9 nor- 
mal stars. 

The BCES algorithm has been applied to the sample of Am plus A3-A9 
normal stars in an iterative procedure, taking into account all the errors in the 
input parameters. The obtained calibration using Hipparcos data is (Domingo, 
1998): 
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M v = (12.16 ± 0.03) - (3.5 ± 0.7)/3 - (8.3 ± 0.3)<5ci + 
+ (6.1 ± 1.0)£mi + (7.1 ± 1.1) • 10" 6 (Fsini) 2 

The unweighted residuals obtained in M v are 0.25 mag for A3-A9 V normal stars 
and 0.36 mag for Am stars, smaller in all cases than for the existing photometric 
calibrations in this spectral region. 
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